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^ ' Abstract 

pL ! 

' , We consider interactions of fcrmions with the domain wall bubbles produced during 

^^ ' a first order phase transition. A new exact solution of the Dirac equations is obtained 

for a wall profile incorporating a position dependent CP violating phase. The reflection 

^ ■ coefhcients are computed, a resonance effect is uncovered for rapidly varying phases. 

T^ I This resonance effect happens when the energy and mass of the incident particles are 

CO ■ E/ra = A6'/2. Where A0 is the phase variation across the wall width. We calculate the 

vJn \ chiral charge flux through the wall surface and the corresponding baryon asymmetry 

of the Universe, ft agrees in sign and magnitude with the observed baryonic excess 

Pb/ s « 10~^° for a large range of parameters and CP violation. As a function of A6', 

^\ • the ratio pb/s reach a maximum for A0 « 2 — 47r and m « rritop- 

PACS: 11.27.+d, 03.65.-w, 02.30.Hq, 02.30. Gp, 11.30.Fs, 98.80.Cq 



1 Introduction 

During the last few years, a considerable amount of work has been dedicated to the possi- 
bility of generating a sizeable baryon asymmetry on the electro-weak phase transition (see 
for example [||, |2|, y , ^ ) . For an excess of baryons to develop in an Universe which initially 
has zero baryon number, the already well known following conditions, first enunciated by 
Sakharov, must be met: 1) Some interaction of elementary particles must violate baryon 
number. 2) C and CP must be violated in order that there is not a perfect equality be- 
tween rates of AB ^ processes, since otherwise no asymmetry could evolve from an 
initially symmetric state. 3) A departure from thermal equilibrium must play an essential 
role, since otherwise CPT would assure compensation between processes increasing and 
decreasing the baryon number. Remarkably, the standard model of weak interactions may 
provide all these conditions. In particular the third one can be met if the weak phase 
transition is at least weakly first order. 

In a first order phase transition, the conversion from one phase to the other occurs 
trough nucleation. This happens when the system is either supercooled or superheated. 
Bubbles of the "true" phase (with an expectation value of some Higgs field v ^ 0) expand 
rapidly absorbing the region of the "false" phase (f = 0). At the bubble surface there is 
a region or "wall", in principle of microscopic dimensions, which separates the phases. 

The speed regime of the expanding bubble walls is poorly understood in the theory. 
After an acceleration period, the final stationary speed of the expanding bubble walls 
could be anywhere in the range 0.1 — 0.9 c P]. Lattice and numerical simulations indicate 
that Vmaii ~ 0.5 can be reached in circumstances of high latent heat and high surface 
tension which could be achieved in the minimal super symmetric standard model (MSSM) 
for realistic Higgs masses MH<^M]y. For having even higher wall speed v^aii ~ 0.9, it 
is necessary some fine tuning (a extremely low value for the friction parameter, see [^). 
At later stages of the bubble evolution an oscillatory regime with successive periods of 
expansion and contraction is naturally expected ^. 

Particles in the "false" (higher temperature) phase are reflected off the advancing 
bubble walls, while most particles in the low temperature phase are unable to catch up 
with the receding walls, and cannot equilibrate across the phase boundary. Thus one has 
a departure from equilibrium and a baryon asymmetry can be generated. 

In the physical conditions of the early Universe the fermions moving through the bubble 
wall will interact also with the particles in the surrounding plasma, a full transport problem 
must be considered. A useful simplifying assumption is to decompose the process into two 
steps, one describing the production of the CP asymmetry on the transmission/reflection 
coefficients when the quarks/antiquarks are scattered on the wall, the second describing the 
transport and the eventual transformation of the CP asymmetry into a baryon asymmetry 
via the baryon number anomaly. 

Assuming that the scattering from the wall is little affected by diffusion corrections, 
the effects of the surrounding plasma can be partially incorporated by introducing a Higgs 
field effective potential which takes into account finite temperature corrections to the tree- 
level potential. The structure and profile of the wall depends on this effective potential 



on a complicated way. Fermions passing through the domain wall acquire a mass which 
is proportional to the vacuum expectation value (VEV) of the Higgs field, which is deter- 
mined from the equations of motion of the finite temperature effective action of the bubble. 
The problem of computation of transmission coefficients reduces to the solution of a Dirac 
equation with a space dependent mass term. Exact solutions has been obtained only for 
two simple cases: for a wall profile approximated by a step function |Q, |^] and for an 
average Higgs field profile of the "kink ansatz" type [^(with z, a coordinate perpendicular 
to the wall): 

(f>iz) oc (l + tanh (^^^ (1) 

where the width of the wall is given essentially by the inverse of the Higgs mass: 



Mh 



(2) 



In this ansatz there is not variation of the complex phase of the Higgs field through the 
wall. It has been examined in [^ the effect of introducing a small CP-violating imaginary 
part in Eq.O as a perturbative effect. Numerical solutions have been obtained for some 
more complicated profiles which incorporate locally dependent complex phases and which 
are considered to be " reasonable" enough although not neccesarily solution of any equation 
of motion ( for example in [^ ) . 

The object of this work is the computation of transmission coefficients of a fermion 
in a certain fixed Higgs field background with an arbitrary, possibly high, complex phase. 
With them we will be able to estimate the baryon asymmetry of the universe. 

Consider a langrangian for the fermion ^ in the background of a complex scalar with 
a Yukawa coupling. For the scalar to have a nonzero complex vacuum expectation value, 
we can assume it is a part of an extended Higgs sector such as MSSM or a more general 
two Higgs-doublets model. The fermionic part of the lagrangian important to us is 

L^ = "^ii jd^L + "^Ri jd^R - h^L'^RCp - h*^R^L4>* (3) 

with Yukawa coupling constant h. ^ acquires a position dependent mass, which, in a 
appropriate reference frame is a function only of one coordinate , m{z) = h < (p{z) >. 
From L^ the equations of motion are obtained 



i jd^L - m*^R = 
i jd'^R - m^L = (4) 

It is possible to absorb a constant phase in m(z) with a redefinition of one of the chiral 
fields. 

One needs only the plane wave solution of the Dirac equation for particles moving along 
the normal to the wall profile (considered microscopically flat). For any other incoming 
direction, the problem can be reduced to the latter performing a Lorentz boost. Reordering 



the spinorial components, the Dh'ac operator given by Eqs.@ can be factorized into 2x2 
blocks. For solutions with positive energy E, we can make the ansatz: 






■04 
VV'2/ 



(5) 



Where xt, is the projection onto the x-y plane. V'l and V'2 are eigenspinors of the chiral- 
ity operator , 75, for the eigenvalue +1 and V'3>'04 ^^r —1. Going through an intermediate 
Lorentz transformation we obtain finally the two equations 

{id, + Q{z))iPj = (6) 

{id, + Q*{z))^ljn = (7) 



Qi^) = {jL ~""i"^);Ei = +\/E^-p't (8) 



With 

n(-y\ — ( 

m{z)* -El 

We will deal in this work with the particular case given by the function 

( \—\ "i-oexpz(— A6'Az) exp — Az if z > , . 

""^^^ ~ I mo if z < ^> 

The constants A0, A are real; A > ; rriQ can be also set to real without loss of generality. 
This function represents a linear phase variation with a global difference of A^ over a 
distance of the order of the wall thickness defined hy 6 = 1/A. In specific models A0 
can be related to the quantity of CP violation in the Higgs sector. The position depen- 
dent complex phase can be considered as an additional source of effective CP-violating 
phenomena. 

The profile of the wall given by m{z) cannot be derived from any classical field equa- 
tions in presence of an effective Higgs potential with up to quartic terms as it is usually 
considered. It can be trivially seen however that such a wall profile corresponds to an 
energy density of the wall per unit area of the form 

£{M) - £:(0) = AmgA02 (10) 

with <f (0) of the same order of magnitude as the energy density of the "kink" ansatz (||). 
Some indication can be obtained from the previous formula. As the nucleation temperature 
Tn oc £ (see Eq.(13) in |jlO|) and according to Table (1) in |p which relates nucleation rates 
and wall speed, a high value for A9 would favor generically a lower value for the wall speed. 
In this work we will see (Fig.(6)(C-D)) that, for small wall speed, a high value for A9 is 
equally acceptable as a lower one in order to obtain a reasonable pb/s ratio. 



2 Solving the Dirac equation. 



It would be possible to solve Eqs.(^0) with the function given by Eq.(^), reducing them 
to a Bessel-like differential equation. It is possible and advantageous however to use 
perturbation theory to compute the evolution operator of the system. The summation to 
all orders of the perturbation expansion is possible thanks to the special form of Q{z). 
There are two main advantages in doing so: the first one is that the procedure is easily 
generalizable to any number of dimensions (for example to incorporate mixing between 
generations), the second one is that the evolution operator is computed directly and the 
reflection coefficients are easily given in terms of its components. A similar technique has 



been used already in [|TT|, 12 1 to compute the neutrino oscillation probabilities in solar 
matter. 

The evolution operator of the differential system (^) is given by the path-ordered 
integral 

U{z,zo) = Pexpi f dzQ{z) (11) 

Jzo 

In this work we are concerned with an operator Q of the form 

Q = Q^ + V{z) (12) 

With 

k real, a complex in general. The real part of a is greater than zero. It will be set 
JRo" = 1 without loss of generality. 

Formally, it is possible to solve Eq.(0) by successive iterations: 

oo 

Uiz,zo) = U^^\z,zo) + ^ f/(")(z,zo); U^''\z,zo) = exp {iQ^iz - zq)) (14) 

n=l 

U^"^' is the well-known integral 

C/W =i^JdZn... dziU\z, Zn)V{Zn) . . . U\z2, z{)V {z{)U\zi , Z^) (15) 

The domain of integration is defined by 

T = Z > Zn > ■ ■ ■ > Zi > Zq. 



Following the same arguments as in |11] one can show that it is enough to compute 



U in the z ^ oo limit. The evolution for finite z can be deduced from the expression for 
this limit. Trough elementary manipulations of Eq.(|l5|), we get the elements of U^'^' in a 
basis of eigenvectors of Q^: 

< b I ;/(") I a >= eexY>{i{Qlz - Q°zo))x 

X Y^ d''rexp{iZnWbki + ■■■ + iziw^n-i)a)Vbkiizn)---V^n~i)aizi) (16) 

kl,...,k{n-^) 



With Wkik2 = Qki ~ Qk2- Qk ^^^ °f ^^^ *^o eigenvalues of Q^. 

Due to the dimensionahty of the problem and the special form for V, the summatory 
in Eq.(16) either is zero or reduces to only one term depending on whether n is odd or 



even and on the states a, b. For diagonal terms (a = b) the product of V . . .V is always 
zero when n is odd. When n is even there is a single surviving term. For non-diagonal 
terms, the situation is reversed: the only one surviving term appears for n odd. This 
single remaining term is always of the alternating form . . . F12V21V12 • • •• So, 

< b I t/W I a >= ." exp(z(QOz - g^o))/^^ x | ^^^ J" | ^ j)"^^!^/, ^ = ^ (17) 

(n) 

All the functions appearing in the integral I^^' are of exponential type, the following 
equality [0] can be applied: 

/•CXD /•X2 

In{wi,... ,Wn) = / .../ dXn-.-dxiexp'^WnXn 
J zn J zn 



-l)''exp(zoEn^n) 



Wi{wi + W2) . ■ ■ {Wi +W2.--+ Wn) 

(valid if 3fJ t^„ < 0, V n ) (18) 

In our case the sums in the denominator of Eq. (|l^) are of the form 

(19) 



wi + W2 + ■.. + Wj = iwbki - L + iwkik2 - L* + ... + iwf,(j^i)kj - L 
= iwbkj + niL + 77-2^* 



'^bk(j) can take only the values {0, ib2£'}. L is o" or cr*. One have for the integers 
ni,n2'- ni + n2 = j ni = n2 or ni = n2 ± 1; they are equal or differ in one unit. 
For the diagonal terms, a = b, n even (writing s = 2Ei — a): 



/("•).' 



,611671 

"■aa 



expzo(n/2)(s + s*) 



s*{s* + s){s* + S + S*)... ((n/2)(s + s*)) 
exp zo(n/2)(s + s*) 

_ ex]3Zo{n/2){s + s*) 

~ {s + s*)"/2 (f )!(s + s*)"/2 utli^Vis + s*) + (/ - 1)] 

_ ex.p zo{n/2){s + s*) 

(s + s*)"(f)'[sV(s + s*)](„,/2) 

For non-diagonal terms, taking a = 1, 6 = 2 to simplify the notation: 



(20) 



/: 



{n),odd _ (-1) exp zo((ra/2)(s + s*) + s*) 



ab 



s*{s* + s){s* + s + s*)... {{n/2){s + s*) + s*) 



-l)expzoin/2){s + s*) 



np2,e.en(j72)(s + S*) U]=1,om{^* + (j " 1)(^ + S*)/2) 

(-l)expzo(ra/2)(s + s*) 

~ (s + s*)("-i)/2 (^) !(s + s*)("+i)/2 Utt^^ ^^[s*/{s + s*) + (/ - 1)] 
(-l)expzo(n/2)(g + s*) 

We have used the Pochammer symbol defined by Eq. (B^). 



So, inserting Eqs.(20,21) in Eq.(17) and taking s + s* = —2, zq = 0: 



(21) 



<l\U\l>-- 






n/2 



1+E 



n=2,ei)en 
fc2 



-2)"(§)![-sV2](„/2) 



;\^4; m![-sV2]M, 






(22) 



< 1 I t/l 2>= 



oo 



n=l,0(id 



(n-l)/2 



(-2)"(^)![-.V2]((n+l)/2) 



JEz — 2«fc y^ 



A;2 



m=l 



-ik 



— oFi l-y;y 



/ (m-l)![-sV2](m) 



and similarly for the other two matrix elements. 

< 2 I C/l 2> 



„—iEz 7-1 

e o-c^i 



2' 4 



<2| [7| 1 > 

The matrix U can be written as 

U{z — > oo, zo) 



iEz 



ik 



s e 



•^■-'-""A^-rT 



ex^iQo{z - ZQ)Ured{zo) 



(23) 

(24) 
(25) 



Ured{zo) 



F G 

Q* p* 



(26) 



writing a = A(l + iAO), k = mo/A, E = Ei/X 

_ -imo/X ll.(Ei-^\i.!li 

1 + {2Ei/\ - M)i''^ I 2 ^ ^ A 2 ^ ' 4A2 



1 , r^_^v-:!i 

2 + ^ A 2 ^ ' 4A2 



0^1 U + (^ - ^)^; 7T^ (27) 



In the case zq 7^ 0, F,G would include a factor exp —2zq in its argument. The special case 
Ei/\ = A0/2 will be important later. In this case F, G become dependent only on rriQ and 
have the explicit simple form: 



° = -'— »-^' U'iA? =""""'VA 



o^Ji;^ =™hf^) (28) 



To obtain U , evolution operator for the Eq.(|^), one must make the change A^ -^ — A6'. 
The matrix become 



with 



G 



t{z^oo,zq) = 


expiQo{z- Zo)Uredizo) 


Ured{zo) = 


(F G\ 
W F*) 


—imo/X 
1 + (2E,/X + Ae 


(3 El Ae ml 
).°^M2 + ^A + 2 ^Sa2 



(29) 



1 ,Ei Ae ml 



See Appendix A for some new formulas for the absolute values of generalized hyper- 
geometric functions which can be obtained from the general properties of U, U. 



Following the same reasoning used in [11|, using the general properties of the evolution 



operator, the evolution for any finite z is given by the matrix 

U{z,zo) ^ Us{z)-'Us{zo) = U~l{z)e''^^^''^'>^Ured{zo) (31) 

In this work, we will make use only of the infinite z limit. 

3 The reflection Coefficients. 

We will follow the same procedure as in [y for defining the reflection coefficient. As it can 
be seen directly from Eqs.(26,29), in the region z — > 00 we have the following asymptotic 
behavior 



Moo), Moo) ~ e'^'^ (32) 

In this region , the eigenspinors '02)V'3 correspond to, vanishing mass, left-moving 
particles of opposite chirality and similarly ipi,ip4 to right moving particles. 

For left-moving particles incident from the symmetric phase two components coexist 
at z — > cxD, the incident particle itself and the reflected one by the domain wall. We define 
asymptotic reflection coefficients R, R as 

'0i(cxd) = Ril;3{oo); ^4(00) = -^^2(00) (33) 

The momentum eigenstates for z < are obtained diagonalizing the constant matrix 
(5(0). The eigenvectors 

correspond to right moving (-|-) and left moving particles (-) with respective eigenvalues 

P± = ±Pi = ^\/Ef -ml. 

Imposing the boundary condition that at z < only a left-moving particle with mo- 
mentum p_ propagates, the wave function at infinite is: 

^(00) = U{z -^ 00, 0)V'- (35) 

The reflection coefficient R is given explicitly by the expression 

and similarly for R. 

We are interested in the asymmetry between the reflection probabilities: A =| i? p 
— I -R p. In the ultrarelativistic limit {Ei/mo » 1), R -^ F* /G. In the limit E/rriQ — > 1 
{t -^ 1), both reflection coefficients | i? |,| i? |^ 1 and the asymmetry A ^ 0. It is 
straightforward to check that in the case we had supposed from the beginning a complex 
constant mass mo, the quantities R, R would remain independent of the phase of ttiq taking 
into account the formula (pq). 

In Fig.(jl]) the asymmetry A is plotted as a function of the dimensionless quantity Ei/rriQ 
and some fixed values of A^. For A0 = vr (plot B) the results obtained here coincide with 
the numerical result obtained in |j^] (Fig. 2). For a realistic case where A ~ Mh ~ M\y 
and rriQ ~ rritap, the value of /i = rriQS = itiq/X ss 2 — 2.5. For this range of fi (continuous 
curves in Fig.(|l])) the asymmetry A is non-neglible for a large E/tuq range. In general A 
is negligible for fermion masses corresponding to /_i<~10^^ which can be translated to the 
condition Tno<^2mbottom for not so heavy Higgs masses. 



As the phase difference A^ increases, the height of the peaks keeps roughly unaltered 
but their position moves: 

(^^Mo)peafc =1 ^^ I /2^- 

For large differences Ei/mQ — /S.6 /2^, R becomes neglible and the asymmetry A ~| R p. 
R is non-neglible only for E/rriQ and A0 small enough, such that the imaginary part of 
the first parameter of the Hyper geometric function in Eqs.(p7D approaches zero. 

For higher values of A0 some resonant effect becomes apparent, the peaks get sharper 
and the asymmetry is only non-negligible around them. This is particularly evident in 
the last plot (D). This effect could lead to an essentially monochromatic asymmetry in 
scenarios with highly oscillatory or random phase differences. For high A0, the asymmetry 
quantity can be approximated by the expression 

A^5{Ei/m^-M/2^)f{^i,M) (37) 

where 6{x) is a delta-like peaked function and /(/i, A/i), the value of A at peak, is easily 
computable from Eqs.(p^), it is of interest only its general behavior: 



The resonant behavior is again evident in Fig. (0) where we plot A as a function of the 
parameter A0 for various // as before. The peaks in the asymmetry function appear for 
values such that A^pgafc = 2£';(5. 

4 Chiral charge flux through the bubble wall and the Baryon 
asymmetry. 

The chiral charge flux in front of the wall in the symmetric phase and in the wall frame is 
given by %. 

1 roo roo 

F = ^:^J^ dpi I Ptdpt{n-Pi,pt) - f{Pi,Pt))A{mo,X,Pi) (39) 

A is the asymptotic asymmetry defined in the previous section and the thermal equi- 
librium fermion flux densities /'*,/'' in the symmetric and broken phases respectively 
(neglecting Fermi blocking factors) are: 

^ exp(7(^- - upi)/T) + 1' ypi^pt 



iPi/E' 



<b\ 



'° = expHg^T^rt/D + l -^'-V^' +"'+'"" 
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The wall velocity in the fluid frame is given by u ; 7 = 1/Vl — "U^; pi,pt are the 
longitudinal and transverse momentum of the fermion with respect to the wall. T is the 
equilibrium temperature. The use of asymptotic expressions for /* and A is questionable, 
the position dependent mass implies the existence of a density gradient even at thermal 
equilibrium. The Formula (^9|) would be valid at far distances from the wall and in absence 
of inter-particle collisions. We can expect that the formula is reasonable valid as long as 
the wall width is small compared to the mean free path in the plasma, so effects inside 
the wall do not play an important role. 

Defining dimensionless momentum and wall width pi/rriQ = x;pt/mQ = y,e = T/X, 
defining also k = fi^/e and integrating with respect pt, we arrive to the expression: 

Where F* is the function defined by the integral: 

F*{^^, K,u)= I dx xAif,, x) I log 1 + exp(-/<l- n)x) j ^^^^ 

1 + exp ( — Kv 1 + x^ — Kux 



The dependence with the CP violating phase A0 is fully contained in the asymmetry 



^(/i, x). The integrand in Eq.(41) is defined positive if u > 0. For negative u (as it could 
happen in late wall contraction phases) its sign changes at a certain xq = l/{2^/u^+ \ u |). 
For large x the integrand multiplying A in Eq.(^) approaches zero exponentially for any 
value of the parameters, the most important contribution to the integral comes from the 
values of A{x) with x relatively small. This implies that in the case F* might become 
negative, its absolute value would be relatively small. It is expected that the resonance 
observed in the previous section tends to be washed out due to this behavior: the resonance 
in A appears at large x, just when the rest of the integrand is exponentially smaller. The 
resonance in A could have a more significant effect in non-equilibrium cases: for flux 
densities f^,f with large tails at large momentum. 

We plot in Fig.(|3|) the chiral flux divided by the cubic temperature given by the Eqs.(40- 
|4l|) as a function of u. F/T'^ is rather independent of u and /S.9 in the logarithmic scale 
of the figure, except for n ^ 1. It is however strongly dependent in ^, it changes a few 
orders of magnitude for ^ varying changing in the range 1 — 10~^. 

For large A0, using the formula (|37|), one obtains the approximate expression for the 
integral F*: 



1 + exp f-K(l - n)v/(A6'/2^)^- 

F* ^ f{fi, A9) log ^-\ ^ v^ (42) 

1 + exp i^-K i^A9/2fi + u^/{A9/2fi)'^ - 1 



For a large value of A9 but keeping the quantity KA9/2fi = 'yA9/2e relatively small, 
F*, and so F, are proportional to A9 

F. = ^/(^,A<,)^ (43) 

11 



This behavior coincides essentially with the proportionality in A0 obtained in [p^ j. 

In the opposite regime, for a large quantity kA0/2//, F* approaches zero exponentially: 

F* ^ 2 fill, Ae)e-''^^/^>' sinh (raA0/2/i) (44) 

We plot F* as a function of A6 in Figs.(^) (given by the the exact Formula (j^). The 
small and large A9 behavior observed in the figure are clearly given by the Eqs.(^^-p3|) 
respectively. It is interesting to compare Fig.(H) with Fig. (pi). In agreement with what 
was expected from the form of the integral (|4l|), the resonances are lost here. However, It 
remains still a non obvious remarkable behavior. F* (and then pfe/s to be defined later) 
is oc A9 in a surprisingly sizeable range, after reaching a maximum for A6 around 2 — 47r 
for practically any value of the rest of the parameters, F* goes quickly to zero for higher 
values. 

4.0.1 The baryon asymmetry 

The baryon asymmetry/entropy ratio of the Universe through the sphaleron process is 
given in the charge transport scenario with two Higgs doublets by |j|]: 

» = i£ii£: (45) 

s b T u s 
Where 



r^ is the sphaleron transition rate per unit of volume in the symmetric phase |14]: 

Tb ^ 3r]atT^ (46) 



with an unknown numerical factor ry ~ 0.1 — 1 [|15|. Note however that a violation 
rate ~ O(a^r^) has been proposed recently [16|. The introduction of an additional 



factor Ouj would not change greatly the final results, or at least its order of magnitude, 
but it is an indication of the global incertitude still existing on F^ . 

T is a typical transport time within which the scattered fermions are captured by 
the wall and in which they can be converted into baryon number |j|] t ~ l/u with 
1 the mean free path of the corresponding particle. In some works r is computed 
numerically Q. We will parametrize it in the form 

T 

D can be interpreted as a diffusion-like constant with a magnitude in the order 
D^l-20 (for example D ?» 4 - 6 in 0, [l|]). 

The entropy density is given by s = 2Tr'^ g^T'^/Ab with g^ ss 100, the effective degrees 
of freedom of the relativistic particles at the electroweak phase transition. 



12 



• F is the charge flux given by Eq.([4(]|) which contains all the dependence on A0 
through the integral F* . 

Inserting all the factors the baryon/entropy ratio is given explicitly by: 

^ « 5.10-9 X {r,D) X 4^ (47) 

s u^e^7^ 

F* is a function only of the adimensional parameters e,/i, n and /S.6. The explicit depen- 
dence with the temperature T has vanished in this formulation. The range of validity 
of Eq.(^) is limited by the simplicity of charge transport model. It should be valid for 
Tbt/T^ « u ([|1|) or 7]D « lO^u^. With r] ^ l,D Ri W, we would have then the 
condition u » 10^^. In the other extreme the formula lack a priori validity for u ^ 1. 
In this case there is not time for the fermion in front of the wall to thermalize and the flux 
densities defined in the previously can not be used. For those cases a detailed Boltzman 
transport treatment should be performed instead. Note however that, according to previ- 
ous arguments taking into account the dependence of A at large momentum, the thermal 
case would act as a lower limit. Another constraint comes from the condition for the 
validity of formula (|39|): the wall width should be smaller that the mean free path of the 
fermion or e = T / \<^D . 

The matter content of the observed universe is given by pb/s ~ 10^^'^ (from nucle- 
osynthesis calculations, see Q and references therein). Figure ^ shows the ratio pb/s. 
For large u, ps/ s follows the behavior of F: decreases and approaches zero for n ^ 1. For 
small u, as long as the equation is still valid, the divergent behavior ~ 1/u^ is apparent. 
According to Fig.(|5|) a value pi,/ s ~ IQ-^*' may be explained by this mechanism for weak 
scale baryogenesis for a wide range of parameters, including a large variation in the rather 
unknown product r]D^ as long as // ~ 0.5 — 2. //is the main parameter needed to be tuned. 

In Figure ^ the same ph/s is plotted as a function of negative wall speed (contracting 
wall case). In agreement with previous considerations ps/ s is of opposite sign in this case 
and relatively smaller in absolute value. The effect is more important for // ~ 1 — 2, A0 
large and medium speeds. Note that the pattern of the p dependence is different here than 
that one appearing in Fig.(|5|). For a detailed knowledge of the final ps/ s an integration 
over the full history of the bubble may be neccesary. 

The dependence with the adimensional thickness of the wall (e = T/X) is shown in 
Fig.(0) ( cf. with Fig. (4) in [Q]). e corresponds approximately to the inverse Higgs mass 
scale and should be independent of T (A ~ T) in most models. The dependence is rather 
modest for thin walls (e < 1). For relatively thicker walls e > 5 — 10 the suppression is 
very important: values as ps/ s « 10^^^ are reached in this region. 

The dependence with the CP violating phase A0 of ph/ s is contained totally in F* 
and appears in Fig.(^. Although interesting qualitatively, it is clear from a numerical 
point of view that the maximum variation , a factor 4 from A0 ~ 7r/2 to A0 ~ 27r, is still 
small in comparison with all the other incertitudes appearing in the present knowledge, 
observational and theoretical, of Pb/s. An important conclusion from Fig.(B) is that, 
although decreasing, pi,/ s still keeps the right order of magnitude even for very large A0. 
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5 Summary and conclusions 

As conclusion, we were able to solve the Dirac equation with a space dependent complex 
mass term which, although not a solution of the equations of motion, reproduces the 
expected variation in module and in complex phase of the average Higgs field across 
the wall. From the solution of the Dirac equation the particle-antiparticle transmission 
asymmetry is computed. The analytical results presented here confirm previous numerical 
computations for small A^ and predict an unexpected behavior for highly oscillatory phase 
fields. 

We have computed the chiral flux and the baryon asymmetry. For a large set of 
parameters this is compatible with the "experimental" value pb/ s ~ lO^^'^. This is true 
even for extremely large values of CP violation IS.6. The resulting baryon asymmetry for 
contracting walls has been estimated for the flrst time to our knowledge. 

The main uncertainties come from the rate of anomalous baryon number violation at 
high temperature and from the method itself. In view of the dependence on e and A^, 
It is clear that the plausibility of this model depends only on the existence of a massive 
fermion (as the Top quark) and is rather independent of the values of the Higgs mass and 
the critical temperature of the weak transition as long as they keep in a reasonable range. 

A more detailed microscopic treatment based in the Boltzman equation is clearly ad- 
visable. The solution of the Dirac equation obtained here is particularly suitable for this 
purpose. Although the mass term which has been used is somehow artificial, is a suit- 
able toy model to study large CP violation environments. This solution can be used as a 
starting point to obtain numerical correction for more realistic mass fiuctuations. 

The method used here can be used with little modifications to solve the Dirac equation 
for a purely local dependent phase mass term (A -^ 0,A9X ^ in Eq.(^) or equivalently 
a purely imaginary in Eq.([T^)). In order to circumvent convergence problems this must 
be done taking the limit JRA ^ in the finite propagation time Eq. (|3l|) . 



A Appendix: some old and new formulas about Hypergeo- 
metric Functions 



The generalized Hypergeometric function ||18[ is defined by 

pFg{ai, ...,ap,h,...,bg;z) = }^ — — — — (48) 

^ (Olj(n) ...(6g)(n) n! 

where the Pochammer symbol is 

(z)(„)=r(n + z)/r(z) (49) 

in particular 

o 
o Fi (h- ^\ = S 

to (^)(n) nl 



CO ^ fl 

o^i(^;^) = E77;^^ (50) 



n= 
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The derivative of o-^i is again a o-^i function: 

d oFi{r,z) 1 77 n I ^ fK^\ 

= -oFi{l+r,z) (51) 

az 7 

The function qFi is related to the Bessel Functions by the formula 

For a Q as given by Eq.([T^) is traceless, det Q{z,zq) = det Q(0,0) = 1, or | F p — | 
G p= 1. The formulas ( |27| ) are vahd for any E, A6 real, k complex, we obtain: 



oFi ( i + e^■, x'^ f "^/^^^ I o^i (^ + «; ^' ) \'= 1 (53) 



to be compared with the similar formulas obtained in |11] for the absolute values of the 
generalized hypergeometric functions nFn- 

In fact, with little changes, one could compute U for any general matrix V with diagonal 
terms equal to zero and non-diagonal terms of general exponential type not neccesarily 
equal. In the particular case of V hermitic: V12 = V21 = k, Ured is unitary and of the form 

Uredizo) = f _Q* p* ] (54) 

with F,G given by Eq.(||). 

By the unitarity of the matrix (|53), |Fp + |Gp=l. And we get the formula: 



J_i+,,(^) P + I Jl^Jx) |2= "^7:"^ (55) 



2 cosh vre 

2 ' "" ■ 2 ' "" ' ■ TTX 

For e = this formula reduces to the trigonometric case involving the well known Bessel 

functions of order 1/2: J-i/2{x) = ^/2J{^^x) cos{x); Ji/2{x) = \/2/(7rx) sin(x). 
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Figure 1: The asymmetry A as a function of Ei/ruQ. Continues line: fJ, = 2, dashed hues: 
respectively fi=l, 1/2, 1/10. From A to D: AO = 7r/2, vr, 2tt, dn. 
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energy is choosen as Ei/rriQ = 2. 
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Figure 3: Chiral charge flux (Eq.([40D) as a function of the wall velocity u, for several 
A^ and fi (continuous and dashed lines as in previous figures). In all the figures the 
adimensional wall width e = 1. 
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rjD = 1) as a function of the wall velocity u (see explanation 
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Figure 7: Dependence with the thickness of the wall e = T/\ . For high and low values of 
the wall velocity and A^. Continuos line: fi = 2, dashed line fi = 1. 
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